3D polarization-interference holographic histology for wavelet-based differentiation of the polycrystalline component of biological tissues with different necrotic states. Forensic applications

Abstract. Significance The interference-holographic method of phase scanning of fields of scattered laser radiation is proposed. The effectiveness of this method for the selection of variously dispersed components is demonstrated. This method made it possible to obtain polarization maps of biological tissues at a high level of depolarized background. The scale-selective analysis of such maps was used to determine necrotic changes in the optically anisotropic architectonics of biological tissues. Objective Development and experimental approbation of layered phase polarimetry of repeatedly scattered fields in diffuse layers of biological tissues. Application of scale-selective processing of the found coordinate distributions of polarization states in various phase sections of object fields. Determination of criteria (markers) for histological differential diagnosis of the causes of necrotic changes in optical anisotropy of biological tissues. Approach We used a synthesis of three instrumental and analytical methods. Polarization-interference registration of laser radiation scattered by a sample of biological tissue. Digital holographic reconstruction and layered phase scanning of distributions of complex amplitudes of the object field. Analytical determination of polarization maps of various phase cross-sections of repeatedly scattered radiation. Application of wavelet analysis of the distributions of polarization states in the phase plane of a single scattered component of an object field. Determination of criteria (markers) for differential diagnosis of necrotic changes in biological tissues with different morphological structure. Two cases are considered. The first case is the myocardium of those who died as a result of coronary heart disease and acute coronary insufficiency. The second case is lung tissue samples of deceased with bronchial asthma and fibrosis. Results A method of polarization-interference mapping of diffuse object fields of biological tissues has been developed and experimentally implemented. With the help of digital holographic reconstruction of the distributions of complex amplitudes, polarization maps in various phase sections of a diffuse object field are found. The wavelet analysis of azimuth and ellipticity distributions of polarization in the phase plane of a single scattered component of laser radiation is used. Scenarios for changing the amplitude of the wavelet coefficients for different scales of the scanning salt-like MHAT function are determined. Statistical moments of the first to fourth orders are determined for the distributions of the amplitudes of the wavelet coefficients of the azimuth maps and the ellipticity of polarization. As a result, diagnostic markers of necrotic changes in the myocardium and lung tissue were determined. The statistical criteria found are the basis for determining the accuracy of their differential diagnosis of various necrotic states of biological tissues. Conclusions Necrotic changes caused by “coronary artery disease–acute coronary insufficiency” and “asthma–pulmonary fibrosis” were demonstrated by the method of wavelet differentiation with polarization interference with excellent accuracy.

The fundamental results obtained within the framework of MMM of biological tissues have wide diagnostic application in various fields of medicine. 37The possibility of obtaining quantitative optical indicators to characterize the evolution of gastric tissue from a healthy state through inflammation to cancer has been demonstrated. 14The criteria for the Mueller-matrix diagnosis of prostate cancer, 38 bowel, 39 and cervix 40 have been determined.A method of differentiation of postmortem traumatic myocardial changes has been developed. 35Polarimetric criteria for determining the time of death have been found. 41he analysis of the literature data has shown that the further successful development of MMM is hindered by two main, as yet unresolved problems.[12][13][14][15][16]42,43 The presence of a high level of depolarized background reduces the depth of modulation (contrast) MMI. 43,44As a result, the sensitivity and accuracy of the polarization diagnosis of pathological conditions of biological tissues decreases. 44,45Studies in this area have revealed the dependence of the depolarization value on the parameters of the optical anisotropy of diffuse biological layers. 46,47To eliminate the effect of depolarization, a model of polar Mueller-matrix decomposition of the biological layer into various components: "polarizerattenuatordepolarizer" was used here.This algorithmic approach made it possible to diagnose cancerous tissues. 48Another promising direction in eliminating the depolarized background may be the polarization-interference layer-by-layer phase selection of scattered laser fields.This approach is used to diagnose diffuse samples of benign and malignant prostate tumors. 49he second problem is that all the data of the image MMM are integrally averaged.11][12][13][14][15][16][17][25][26][27][28][29][30][31][50][51][52] This circumstance limits the functionality of the MMM.3][34][35][36] Scale-selective wavelet analysis [53][54][55][56][57][58] of polarization maps and MMI of optically thin layers of biological preparations of tissues and liquids is used here.][34] Polarization differentiation of formed tumors is more effective in analyzing transformations of the amplitudes of wavelet coefficients for large-scale structures of biological crystals. 35,36hus, for the further development of the imitative MMM of diffuse biological tissues, it is relevant to combine the considered "depolarization" and "wavelet" methods into a single polarimetric technology.
In our work, we propose a new method for implementing this task.Our method is based on the principles of phase scanning of repeatedly scattered object fields of biological tissues.For this purpose, polarization-interference registration of the laser object field is carried out.According to the obtained interferograms, the digital holographic method reconstructs the distributions of complex amplitudes of such a field.The phase scanning of such a field ensures the selection of variously scattered components.In the phase plane of single scattering, the polarization maps of the optically anisotropic architectonics of biological tissue are algorithmically determined.For the found azimuth and ellipticity distributions of polarization, the wavelet decomposition algorithm is used.4][55][56][57][58] At each scale (b) of the MHAT function, the linear dependences (a) of the amplitudes of the wavelet coefficients of the azimuth maps and the ellipticity of the polarization are determined.Based on the obtained dependencies, the central statistical moments of the first and second orders are calculated, which characterize the mean and variance of fluctuations of different-scale amplitudes of the wavelet coefficients.Thus, the most sensitive diagnostic markers of structural changes in the optically anisotropic architectonics of diffuse biological tissue are determined.Further, within the framework of evidence-based medicine algorithms, the operational characteristics of the diagnostic power of the method are calculated.
Our work is aimed at the fundamental development and experimental approbation of this method of polarization-interference polarimetry of repeatedly scattered fields.As objects, we considered optically thick depolarizing samples of native histological sections of myocardium and lung tissue with different optically anisotropy architectonics.
The applied aspect of the work consisted of determination of criteria (markers) for differential diagnosis of necrotic changes in biological tissues with different morphological structure.Two cases are considered.The first case is the myocardium of those who died as a result of coronary heart disease (CHD) and acute coronary insufficiency (ACI).The second case is lung tissue samples of deceased with bronchial asthma (BA) and fibrosis.

Stokes polarimetry of the object field
We chose a right-circularly (⊗) polarized laser beam as the radiation illuminating biological tissues.This condition is necessary when measuring a series of samples of biological tissues.For other states of polarization, the result of its object transformation will be azimuthally dependent on the rotation of the sample relative to the direction of irradiation. 44,45he Stokes vector of such beam has the following form: 4,5,[7][8][9][10][11][12] E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 1 ; 1 1 7 ; 1 3 5 Ushenko et al.: 3D polarization-interference holographic histology. . .

Single scattering
For the case of single scattering (j ¼ 1), the polarization properties of the linear birefringence protein fibrilla in the point with coordinate (r) correspond to the Mueller matrix operator 6,19-25,50-52 ðfWgÞðrÞ j¼1 : E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 2 ; 1 1 4 ; 6 8 1 Here, ω lk is the elements of the Mueller matrix with column (l ¼ 1; 2; 3; 4) and row (k ¼ 1; 2; 3; 4) indexes; ρ is the orientation of the optical axis of the birefringent fibril; δ ¼ 2π λ Δnd is the phase shift; Δn is the linear birefringence index LB; d is the geometrical size; and λ is the wavelength.
The process of single transformation in local point (r) of the probing beam S 0 ð⊗Þ is described by the following matrix equation: e m p : i n t r a l i n k -; e 0 0 3 ; 1 1 4 ; 4 6 6 ðS Ã ðrÞÞ j¼1 ¼ fWgðrÞ j¼1 S 0 ð⊗Þ: Here, ðS Ã ðrÞÞ j¼1 is the Stokes vector of a single scattered component of an object field at a point ðrÞ: E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 4 ; 1 1 4 ; 4 1 7 where αðrÞ is the azimuth and βðrÞ is the ellipticity of polarization.
Thus, at each point r repeatedly (j ≥ 1; 2; 3; : : : p − 1; p) of the scattered field, the azimuth and ellipticity of polarization are averaged to certain values α p ðrÞ and β p ðrÞ: E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 0 ; 1 1 7 ; 5 5 8 α j¼p ðr; ρÞ ¼ 0.5 arctan E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 1 ; 1 1 7 ; 5 0 3 As a result, a polarization-inhomogeneous component of the diffuse field will be formed with a different distribution of azimuth values A j¼p ðα j¼p ðR; ρÞÞ and ellipticity B j¼p ðβ j¼p ðR; δÞÞ polarization.Thus, a polarization-structural map of the diffuse layer of birefringent biological tissue is formed E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 2 ; 1 1 7 ; 4 The conducted consideration [ratios (1)- (12)] describes the process of direct formation by optically anisotropic architectonics of the "object component" PðR; α; βÞ of a polarization inhomogeneous field In parallel with the direct acts of interaction of laser radiation with birefringent biological crystals [ratios (1)-( 12)], secondary interference of coherent scattered waves occurs.The result of this process is the amplitude addition of variously polarized (ratio ( 13)) partial coherent waves.

Amplitude consideration
8][29][30][31] Based on this, the expressions ( 6) and ( 7) obtained earlier for the polarization parameters can be rewritten in the following form: E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 5 ; 1 1 7 ; 1 1 8 Here, i is an imaginary unit and * is a complex conjugate quantity.

Interference interaction
For the orthogonal amplitude components U j¼1 x and U j¼1 y , the following interference equations can be written: E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 6 ; 1 1 4 ; 6 9 1 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 7 ; 1 1 4 ; 6 4 4 where jU j¼1 x1 j; jU j¼1 y1 j are the modules of complex amplitudes; δ x12 и δ y12 are the values of phase shifts between ðU j¼1 x1 ; U j¼1 x2 Þ and ðU j¼1 y1 ; U j¼1 y2 Þ.For the process of forming the resulting values of the orthogonal components of the amplitudes U x and U y as a result of the j ¼ p interaction, we can write The interference addition of two phase-shifted φ p xy orthogonal components U j¼p x and U j¼p y forms an elliptically polarized wave 6 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 9 ; 1 1 4 ; 4 9 9 with next interference means (I) of azimuth α j¼p ðIÞ and ellipticity E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 1 ; 1 1 4 ; 4 1 2 Thus, a secondary polarization-inhomogeneous component [ratios (20) and ( 21)] of a diffuse object field with probabilistic distributions of azimuth and ellipticity values is interferentially formed ; t e m p : i n t r a l i n k -; e 0 2 2 ; 1 1 4 ; 3 6 9 PðIÞ j¼p ðR; αðIÞ j¼p ; βðIÞ j¼p Þ ⇔ AðIÞ j¼p ðαðIÞ j¼p ðR; ρÞÞ BðIÞ j¼p ðβðIÞ j¼p ðR; δÞÞ :

Resulting field
So, the polarization structure of the laser field diffusely scattered by the tissue layer can be represented as a superposition of the following components:

Wavelet analysis of polarizing maps
The wavelet decomposition provides the mathematical possibility of large-scale selective analysis of a polarization-inhomogeneous field [relation (23)].A salt-like MHAT function with a variable scale (b) of the coordinate scanning window (a) is used as a wavelet.][34][35][36] The continuous wavelet transforms of the function Φðα; βÞ is defined by the following equation: [53][54][55][56][57][58] E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 4 ; 1 1 4 ; 1 2 1 where a is a scale parameter, b is a spatial coordinate, and Ψ is a soliton-like function (wavelet) constructed on the basis of derivatives of the Gaussian function.
The wavelet relations ( 24) and ( 25) for polarization maps of azimuth and ellipticity [ratio (23)] can be written as the following expressions: E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 6 ; 1 1 7 ; 6 6 3

<
: The analysis of expression (26) shows that the wavelet analysis of polarization maps is a superposition of a large-scale selective evaluation of the components of the object field with different light scattering multiplicity.In the following, we focus on the search for the possibilities of optical separation of the components of the object field with a low scattering multiplicity or scattered once.

Experimental Setup and Measurement Methodology
The methodology for polarization-interference measurement of the distributions (m; n-CCD pixels quantity) of the Stokes vector parameters [polarization maps αðm; nÞ and βðm; nÞ] is presented in Refs.27, 28, 35, 36, 38, 44, 45, 49.However, detailed information is not provided in this work.For a better understanding of the further discussion, we provide a brief overview of the three-dimensional (3D) digital holographic scanning method.
A generalization of the polarization interferometry scheme 27,28 is the Stokes-polarimetric mapping scheme on the base of Mach-Zehnder interferometer, which is shown in Fig. 1.
The "object" beam with the help of a rotating mirror 5 is directed through the polarizing filter 6 -7 (manufacturer: Achromatic True Zero-Order Waveplate and manufacturer: B + W Kaesemann XS-Pro Polarizer MRC Nano) in the direction of the biological layer 8 sample.The polarization-inhomogeneous image of biological tissue histological Sec. 8 is projected by the strain-free objective 12 (manufacturer: Nikon CFI Achromat P, focal length: 30 mm, numerical aperture: 0.1, magnification: 4×) into the digital camera 14 [The Imaging Source DMK 41AU02.AS, monochrome 1/2 "CCD, Sony ICX205AL (progressive scan) resolution: 1280 × 960; size of the photosensitive area: 7600 × 6200 μm; sensitivity: 0.05 lx; dynamic range: 8 bit, SNR: 9 bit); the photosensitive area of which contains m × n ¼ 1280 × 960 pixels) plane.The "reference" beam is directed by the mirror 4 through the polarization filter 9-10 (manufacturer: Achromatic True Zero-Order Waveplate and manufacturer: B + W Kaesemann XS-Pro Polarizer MRC Nano) into the polarization image plane of biological tissue histological section 8.
As a result, an interference pattern is formed, the coordinate intensity distribution of which is recorded by a digital camera 14 through a polarizer 13.
Before carrying out measurements of biological tissues, the experimental device passed metrological certification with the introduction of model objects ("clean air," "linear polarizer," "phase plates 0.25λ," "0.5λ").As 50 measurements result for each type of object, the polarization ellipticity errors were determined β ¼ 0.0003 rad.
The object field is scanned to such a value of the phase shift δ ⋆ t , starting from which the statistical condition of a single scattering in the volume of biological tissue is realized: E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 3 1 ; 1 1 7 ; 7 1 3 Z i¼1;2;3;4 ðαðδ t ≤ δ ⋆ t ÞÞ ≈ const and Z i¼1;2;3;4 ðβðδ t ≤ δ ⋆ t ÞÞ ≈ const: (31)

Implement the wavelet transform of polarization maps [ratio (23)]
E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 3 2 ; 1 1 7 ; 6 7 5 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 3 3 ; 1 1 7 ; 6 2 5 9. For various a MHAT scales, the mean Z 1 and variance Z 2 of coordinate dependencies are calculated: ; t e m p : i n t r a l i n k -; e 0 3 4 ; 1 1 7 ; 5 7 9 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 3 5 ; 1 1 7 ; 5 2 6

Objects of Investigations
To implement the complex study of the myocardium histological sections samples: • CHD: group 1, 12 samples; • ACI: group 2, 12 samples.
The second two groups consisted of lung tissue histological sections of those who died from myocardial infarction and with the following concomitant pathologies:   The analysis of images of lung tissue with BA shows a weakly structured parenchymal morphological structure [Fig.2(c)].The image of the tissue with PF shows the presence of a formed fibrillar connective tissue component with an average fiber size of 5 to 10 μm [Fig.2(d)].
This choice of objects is related to both the fundamental and applied components of our research.

Fundamental component
From the physical point of view, the selected tissues have different optically anisotropic architectonics.
The myocardium is characterized by spatially structured linearly birefringent [ratio (2)] networks of myosin fibrils. 5,22,23,35he lung tissue is predominantly parenchymal with a slight linear birefringence of the collagen fibrillar networks of the connective tissue collagen component. 5,29ases of CHD and ACI lead to various necrotic changes of optical anisotropy in myosin fibers and their spatially structured networks.CHD-myosin fibers thin out and the spatial order of fibrillar networks decreases [Fig.2(a)].As a result, the magnitude of the structural anisotropy decreases.ACI-myosin fibers are broken in some areas with constant special ordering of the music network [Fig.2(b)].Therefore, the birefringence level remains commensurate with the same anisotropy parameter of healthy tissue.BA does not significantly change the architectonics of collagen fibrillar networks of the connective tissue component [Fig.2(c)].PF is accompanied by a significant growth of collagen fibers and an increase in birefringence [Fig.2(d)].
Therefore, a comparative physical analysis of the results obtained will allow us to determine the capabilities of our method in detecting changes in the optically anisotropic architectonics of biological tissue samples with different morphologies.

Applied component
An important and not fully solved forensic task by light microscopy methods is to determine the natural (CHD) or violent (ACI) cause of myocardial death.The main problem for histological diagnosis is the presence of a high level of diffuse background in the images of histological sections of the myocardium.
The optical and geometric parameters of the diffuse histological section samples are presented in Table 1.
The extinction coefficient ðτ; cm −1 Þ of the samples of biological tissues was measured according to the standard method of photometry of the attenuation 59 of the intensity of the illuminating beam by the sample using an integral light-scattering sphere. 6025][50][51][52] Histological sections were prepared using the conventional technique on a microtome with rapid freezing. 29ble 1 Optical and geometric parameters of histological section tissues samples of both types.

Information Analysis
For the myocardium, our main applied task was to determine the possibility of detecting ACI cases at a high level of depolarized background.In this sense, the samples of histological sections of the myocardium of those who died from CHD formed control group 1. Accordingly, group 2 of myocardial samples of those who died from ACI was experimental.
The situation was different with lung tissue differentiation.The diagnostic purpose was to determine FP.Therefore, histological sections of lung tissue with BA formed the control group 3, and tissue samples with fibrosis formed the experimental group 4.
Information analysis of the results of polarization-interference phase scanning uses a number of operational characteristics of evidence-based medicine: 61 • Sensitivity (Se) is the proportion of correct positive results (A) of the diagnostic method among all samples from experimental group 2, group 4 (N): ; t e m p : i n t r a l i n k -; e 0 3 6 ; 1 1 7 ; 5 8 4 Se ¼ ðA=NÞ100%: • Specificity (Sp) is the proportion of correct negative results (B) of the method among control group 1, group 3 (H): ; t e m p : i n t r a l i n k -; e 0 3 7 ; 1 1 7 ; 5 3 4 Sp ¼ ðB=HÞ100%: • Accuracy (Ac): proportion of corrPect test results ðA þ BÞ among all samples ðN þ HÞ E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 3 8 ; 1 1 7 ; 4 9 6 If ðA þ BÞ ¼ ðN þ HÞ, then Ac is called the balanced accuracy.
In our work, the following scale for evaluating the diagnostic accuracy is used (Table 2).

Results and Discussion
By the method of phase scanning [ratios ( 27)-( 29)], we have determined the phase plane of a single scattering [ratio (31)].All further discussion refers to the experimental results obtained in this plane δ ⋆ t ≤ π=8.

Wavelet Differentiation of Azimuth Polarization Map of the Myocardium
Figure 3 shows the maps and histograms of distributions of random values of the azimuth of polarization of object fields of the myocardium from group 1 and group 2.
The analysis of the obtained data revealed the presence of a coordinate-heterogeneous [Figs.3(a   Figure 4 shows the results of wavelet analysis of the polarization azimuth αðδ ⋆ t ≤ π=8; m; nÞ maps of the myocardium histological sections microscopic images.
This drawing consists of two parts.On small scales of polarization maps, such wavelet correlations can be traced for fine-structured elements of myocardial fibrillar architectonics.For large scales, variations in the amplitude of the wavelet coefficients are associated with the degree of ordering of large-scale fibers of the fibrillar network.
In the case of death as a result of ACI, the fibrillar myosin network is more ordered and broadly structured in comparison with the architectonics of the myocardium of those who died from CHD [Figs.2(a Quantitatively, we determined the differences between the polarization maps of the myocardial object fields in two ways. The first is a direct calculation method [ratio (30)] statistical moments that characterize the distributions (Table 3).
From the analysis of the data obtained, it can be seen that the maximum accuracy of the differential diagnosis of CHD and ACI practically achieved only a satisfactory level AcðZ 3;4 ¼ 79.2%Þ.
The second is a method for estimating fluctuations in the amplitude of the wavelet coefficients by calculating [ratios (31) and (32)] statistical parameters [mean Z 1 ðW a;b Þ and variation Z 2 ðW a;b Þ] for all scales of the MHAT function (Fig. 5).Here, based on the obtained statistical analysis data, the scales for which the differences were determined.The differences between Z 1 ðW a;b Þ and Z 2 ðW a;b Þ were maximal (Table 4).From the obtained data, it was revealed individual modulation of the amplitude values W b ða min ¼ 15Þ and W b ða max ¼ 55Þ.
For the case of CHD, the maximum amplitude and variations in the wavelet value of the Wa, b coefficients occur on large scales of the wavelet function scanning in the phase plane δ ⋆ t ¼ π=8 of the αðδ ⋆ t ¼ π=8; m; nÞ map.For the case of ACI, the maximum amplitude and variations in the wavelet value of the W a;b coefficients occur on small scales of the wavelet function scanning in the phase plane δ ⋆ t ¼ π=8 of the αðδ ⋆ t ¼ π=8; m; nÞ map.Analysis of statistical estimation of fluctuations in the amplitudes of multiscale wavelet coefficients W a;b distributions of the polarization azimuth αðδ ⋆ t ¼ π=8; m; nÞ revealed next ratios As a result, the following levels of accuracy of differential diagnosis of cases were determined CHD and ACI.For small scales, a min ¼ 15-very good AcðZ 1 Þ ¼ 91.7% and excellent AcðZ 2 Þ ¼ 95.8% are equal.For large scales, a max ¼ 50 is a good level of AcðZ 1;2 Þ ¼ 87.5%.

Wavelet Differentiation of the Ellipticity Polarization Maps of Myocardium
Figure 6 shows maps and histograms of distributions of random values of the magnitude of ellipticity of polarization of object fields of histological sections of the myocardium from group 1 and group 2.
As for the polarization azimuth maps (Fig. 3), the coordinate distributions of ellipticity turned out to be coordinate-inhomogeneous [Figs.6(a It can be seen that the set of polarization domains of ellipticity maps [ratios (7), ( 8), ( 11), (12)] is individual for the object fields of myocardial samples from different groups.At the same time, statistical distributions of random values of the magnitude of the ellipticity of polarization do not carry information about the large-scale (topological) structure of necrotic changes in the architectonics of birefringent fibrillar networks of the myocardium.
Table 5 presents the results of a statistical analysis of the distributions of the ellipticity of polarization of myocardial samples from the groups CHD and ACI.
It was found that the maximum accuracy of differential diagnostics of CHD and ACI exceed a satisfactory level AcðZ 4 ¼ 83.3%Þ.
Figure 7 shows the results of the application wavelet analysis technique Wða; bÞ for maps of the ellipticity βðδ ⋆ t ¼ π=8; m; nÞ.We obtained fluctuating surfaces of the amplitudes of the wavelet coefficients At all scales of polarization maps of ellipticity, the wavelet correlations are determined by the level of structural anisotropy or birefringence.The value of this optical anisotropy parameter is related to the degree of spatial ordering of the fibers of the fibrillar network.
Therefore, in case of death, as a result of ACI, more amplitude values are formed and the depth of modulation of their changes is greater [Figs.4(a 4).
High levels of differential diagnosis accuracy have been established for small scales a min ¼ 15 is an excellent level of AcðZ 1 Þ ¼ 95.8% and AcðZ 2 Þ ¼ 100%.For large scales a max ¼ 50 is a very good level of AcðZ 1;2 Þ ¼ 91.7% (Table 6).

Wavelet Differentiation of the Polarization Azimuth Maps of Lung Tissue
Figure 9 shows maps and histograms of distributions of random values of the azimuth of polarization of object fields of histological sections of parenchymal lung tissue from group 3 and group 4.  As a result, the differences between the values of the statistical moments of the first to fourth orders Z 1;2;3;4 are insignificant and do not exceed 15% to 25% (Table 7).
From the analysis of the data obtained, it can be seen that the maximum accuracy of the differential diagnosis of CHD and ACI does not reach a satisfactory level AcðZ 1;2;3;4 < 80%Þ.In the case of BA, the optical birefringence of the pulmonary parenchyma is insignificant [Fig.2(c)].For the PF situation, the level of structural anisotropy increases due to the proliferation of connective tissue [Fig.2(d)].
Therefore, as a result of the wavelet transform of the polarization maps of the sample from group 3, small amplitude values with a small modulation depth are formed in comparison with the data for the sample from group 4 [Figs.10(a), (b) and (c), (d)].
We quantified such fluctuations in the amplitude of the wavelet coefficients W a;b by calculating statistical parameters (Z 1 ðW a;b Þ) and variation Z 2 ðW a;b Þ for all scales a of the MHAT function Ψð x−b a Þ (Fig. 11).The scales proved to be diagnostically optimal a min ¼ 22 and a max ¼ 43 (Table 8).
High levels of differential diagnosis accuracy have been established.For small scales, a min ¼ 22 is a good and very good level of AcðZ 1 Þ ¼ 90% and AcðZ 2 Þ ¼ 93.3%.For large scales, a max ¼ 43 is a very good AcðZ 1 Þ ¼ 93.3% and excellent AcðZ 2 Þ ¼ 96.7% level.
The obtained results can be physically related to changes in all scales of geometric dimensions of parenchymal connective tissue birefringence fibrils in the lung parenchyma volume in the case of fibrosis.Therefore, at all scales, there is a maximum modulation of the polarization azimuth αðδ ⋆ t ¼ π=8; m; nÞ and, accordingly, variations in the amplitudes of the wavelet coefficients W a;b .For BA, the modulation of the wavelet coefficients is minimal.

Wavelet Differentiation of Ellipticity Polarization Maps of Lung Tissue
Figure 12 shows maps and histograms of distributions of random values of the magnitude of ellipticity of polarization of object fields of histological sections of the lung tissue from group 3 and group 4.
As for the polarization azimuth maps (Fig. 9), the coordinate distributions of the ellipticity of the object field of lung tissue samples turned out to be coordinate-inhomogeneous [Figs.12(a As a result, the differences between the values of statistical moments of the first to fourth Z 1;2;3;4 orders are 25% to 35% (Table 9).
It can be seen that the maximum accuracy of the differential diagnosis of BA and PF corresponds to a satisfactory level AcðZ 1;2;3;4 Þ ¼ 80% to 83.3%.At the same time, the differences in the topological structure of the polarization ellipticity maps are more pronounced for the fields of histological sections of lung tissue from group 3 and group 4.
Figure 13 shows wavelet maps Wða; bÞ of βðδ ⋆ t ¼ π=8; m; nÞ for BA and PF cases.As can be seen (Table 10), the value of the mean and variance of linear dependencies of the wavelet coefficients of ellipticity of maps of lung tissue samples from group 4 (PF) is three to four times greater than similar statistical parameters for the case of BA.Thus, high levels of accuracy in the differential diagnosis of cases of BA and PF have been established.An excellent accuracy level of 96.7% to 100% was obtained on all a scales of MHAT function Ψð x−b a Þ (Fig. 14).

Conclusions
1.The polarization-interference method of mapping and phase selection of diffuse layers of biological tissues with different morphological structure scattered with different multiplicities of polarization-inhomogeneous components of the object field is analytically substantiated and experimentally tested.2. By the method of digital holographic reconstruction with phase scanning of complex amplitude distributions, polarization maps of a single scattered component in the object field of diffuse depolarizing histological sections of fibrillar myocardium and parenchymal lung tissue with the following types of pathology were algorithmically obtained: • myocardium: CHD -ACI; • lung tissue: asthma (BA)fibrosis (PF).3. Using the scanning soliton-shaped MHAT function, a scale-selective wavelet decomposition of a series of polarization maps of azimuth and ellipticity of polarization is implemented.The geometric scales of the structural elements of the polarization maps of a single scattered component in diffuse fields have been determined and physically justified to differentiate pathological changes in the optically anisotropic architectonics of myocardial and lung tissue samples.4. For the obtained scales, within the framework of statistical analysis of linear distributions of the magnitude of the amplitudes of the wavelet coefficients, markers of differential diagnosis are determined CHD-ACI and BA-PF. 5. Diagnostic relationships between first-and fourth orders statistical moments, which characterize the wavelets coefficients distributions of azimuth and ellipticity maps, and diagnostic levels of various pathological conditions differentiation accuracy are revealed.For cases "CHD-ACI" installed very good AcðZ 1;2 ; a min ¼ 15Þ ¼ 91.7% excellent level of AcðZ 1 Þ ¼ 95.8% and AcðZ 2 ; a max ¼ 50Þ ¼ 100%.For cases "BA-PF" installed very good AcðZ 1 ; a min ¼ 22Þ ¼ 92.3% and excellent AcðZ 2 ; a min ¼ 22; a max ¼ 43Þ ¼ 100% levels.

Disclosures
The authors declare no conflicts of interest.

8 U
y;Ω¼0 deg;90 deg is the orthogonal components of complex amplitude; Ã denotes the complex conjugation operation; and ðυ; νÞ are the spatial frequencies.4. The results of the digital Fourier transform (relation (27)) are used to obtain complex amplitudes distributions according to the following algorithms: E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 8 ; 1 1 4 ; 3 7 0 deg ðm; nÞ → jU ⊗ x;Ω¼0 deg jðm; nÞ; (28) E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 2 9 ; 1 1 4 ; 3 4 0 U 90 deg ðm; nÞ → ðjU ⊗ y;Ω¼90 deg j expðiðδ ⊗ x − δ ⊗ y ÞÞÞðm; nÞ:
) and 3(c)] and statistically distributed [Figs.3(b) and 3(d)] structure of polarization azimuth maps for object fields of myocardial samples from both groups.Topologically, azimuth polarization maps are formed by local domains with different geometric dimensions [Figs.3(a) and 3(c)].The distribution of the scales of such polarization domains [ratios (2) to (7)] is individual for the object fields of myocardial samples from different groups.This can be explained by the maximum correlation in the phase plane of a single scattering of the scales of optically anisotropic fibrillar structures of the myocardium and polarization domains.At the same time, the histograms of the distributions [Figs.3(b) and 3(d)] are quite similar in localization of extremes and ranges of variation of random values of the azimuth of polarization.In other words, statistical distributions of polarization parameters do not provide information about the large-scale (topological) structure of necrotic changes in myocardial architectonics.

Figures 4 (
a) and 4(c) show a two-dimensional array of values of the amplitudes of the wavelet coefficients.It is formed as a result of the algorithmic [ratios (24) to (26)] line-by-line scanning of MHAT distributions by a variable-scale function Ψð x−b a Þ.Figures 4(b) and 4(d) show an example of linear distributions of the amplitude of the wavelet coefficients for two scales of the MHAT function.It can be seen from the data obtained that the 2D amplitude distributions of the wavelet coefficients represent W a;b a fluctuating surface [Figs.4(a) and 4(c)].The amplitudes and the oscillation period are distributed differently for different scales a of the salt-like MHAT function Ψð x−b a Þ.As the scale increases ða ↑Þ, the oscillation period of the linear distributions of the wavelet coefficients W b increases.The amplitude modulation depth of the wavelet coefficients W b is individual for different scales [Figs.4(b) and 4(d)].The revealed patterns can be explained by the following considerations.The magnitude of the amplitude of the wavelet coefficient at each scanning point is determined by the degree of mutual correlation of the geometric parameters of the polarization domain [Figs.3(a) and 3(c)] and the scale of the "window" of the MHAT function Ψð x−b a Þ.The greater the cross-correlation, the greater the amplitude of the wavelet coefficient.

Fig. 3
Fig. 3 (a), (b) Coordinate and (c), (d) probabilistic distributions of random values of the azimuth of polarization of object fields of myocardial samples of deceased as a result of (a), (c) ACI and (b), (d) CHD.
) and 2(b)].Therefore, as a result of the wavelet transform of the polarization maps of the myocardial object field for the sample from group 2, more amplitude values and a greater depth of modulation of their changes are formed [Figs.4(a), 4(b) and 4(c), 4(d)].

Fig. 4
Fig. 4 Maps (left column) and multi-scale linear cross-sections (right column) of the polarization azimuth αðδ ⋆ t ¼ π= 8; m; nÞ wavelet coefficients W a;b of myocardium histological sections of those who died from CHD (top row bottom row) and ACI (bottom row).Two-dimensional array of values of the amplitudes of the wavelet coefficients (a), (c) and linear distributions of the amplitude of the wavelet coefficients for two scales of the MHAT function (b), (d).

Fig. 5
Fig. 5 Scale dependences of (a) the mean Z 1 and (b) the variance Z 2 of the wavelet coefficients W a;b of the αðδ ⋆ t ¼ π= 8; m; nÞ.Explanation in the text.
W a;b [Figs.7(a) and 6(c)].The amplitudes and the oscillation period of the linear distributions of the wavelet coefficients increase as the scale a of the MHAT function Ψð x−b a Þ increases.The depth of modulation of the amplitude of the wavelet coefficients of the polarization ellipticity maps (Fig. 6) is individual for different scales [Figs.7(b) and 7(d)].
), 4(b) and 4(c), 4(d)].We quantified such amplitude fluctuations of the wavelet coefficients W a;b by calculating statistical parameters [Z 1 ðW a;b Þ and Z 2 ðW a;b Þ] for all scales a of the MHAT function Ψð x−b a Þ (Fig. 8, Table

Fig. 6
Fig. 6 (a), (c) Coordinate and (b), (d) probabilistic distributions of random values of the magnitude of the ellipticity of polarization of object fields of myocardial samples that died as a result of (a), (b) CHD and (3), (4) ACI.

Figure 10 presents
the results of wavelet analysis of polarization maps αðδ ⋆ t ¼ π=8; m; nÞ for BA and PF cases.The presence of individual fluctuations in the amplitude of the wavelet coefficients at all scales of the salt-like MHAT function was found [Figs.10(a) and 10(c)].As in the case of wavelet decomposition of polarization maps of myocardial samples, the magnitude of the amplitude at each scanning point is determined by the degree of mutual correlation of the geometric parameters of the polarization domain [Figs.9(a) and 9(c)] and the scale of the "window" of the MHAT function.

Fig. 7
Fig. 7 Maps (left column) and multi-scale linear cross-sections (right column) of the polarization ellipticity βðδ ⋆ t ¼ π= 8; m; nÞ wavelet coefficients W a;b of myocardium histological sections of those who died from CHD (top row) and ACI (bottom row).Two-dimensional array of values of the amplitudes of the wavelet coefficients (a), (c) and linear distributions of the amplitude of the wavelet coefficients for two scales of the MHAT function (b), (d).

Fig. 8
Fig. 8 Scale dependences of (a) the mean Z 1 and (b) the variance Z 2 of the wavelet coefficients W a;b of the βðδ ⋆ t ¼ π= 8; m; nÞ.Explanation in the text.

Fig. 9
Fig. 9 (a), (c) Coordinate and (b), (d) probabilistic distributions of random values of the azimuth of polarization of object fields of lung tissue samples with (a), (b) BA and (c), (d) PF.
) and 12(c)].The histograms of the distributions of random values of the ellipticity of the object fields of the samples from group 3 and group 4 are quite similar [Figs.12(b) and 12(d)].

Fig. 10
Fig. 10 Maps (left column) and multi-scale linear cross-sections (right column) of the polarization azimuth αðδ ⋆ t ¼ π= 8; m; nÞ wavelet coefficients W a;b of lung tissue histological sections those who died from BA (top row) and PF (bottom row).Explanation in the text.

Fig. 11
Fig. 11 Scale dependences of (a) the mean Z 1 and (b) the variance Z 2 of the wavelet coefficients distributions for polarization azimuth maps αðδ ⋆ t ¼ π= 8; m; nÞ.Two-dimensional array of values of the amplitudes of the wavelet coefficients (a), (c) and linear distributions of the amplitude of the wavelet coefficients for two scales of the MHAT function (b), (d).

Fig. 12 (
Fig. 12 (a), (c) Coordinate and (b), (d) probabilistic distributions of random values of the magnitude of the ellipticity of polarization of object fields of lung tissue samples of deceased as a result of (a), (b) BA and (c), (d) PF.

Fig. 13
Fig. 13 Maps βðδ ⋆ t ¼ π= 8; m; nÞ (left column) and multi-scale linear cross-sections (right column) of the polarization ellipticity wavelet coefficients W a;b of lung tissue histological sections those who died from BA (top row) and PF (bottom row).Explanation in the text.

Fig. 14
Fig. 14 Scale dependences of (a) the mean Z 1 and (b) the variance Z 2 of wavelet coefficients distributions for βðδ ⋆ t ¼ π= 8; m; nÞ.Explanation in the text.

Table 2
Threshold levels of balanced accuracy.

Table 4
Mean Z 1 and variance Z 2 of the wavelet coefficients W a;b of the αðδ ⋆ t ¼ π= 8; m; nÞ.

Table 8
Mean Z 1 and the variance Z 2 of the wavelet coefficients distributions for polarization azimuth maps αðδ ⋆ t ¼ π= 8; m; nÞ. a min ¼ 22

Table 10
Mean Z 1 and the variance Z 2 of the wavelet coefficients distributions for polarization azimuth maps βðδ ⋆ t ¼ π= 8; m; nÞ. a min ¼ 22